
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



THE LAPLACE-POISSON MIXED EQUATION. 

By K. P. Williams. 

In a recent number of this Journal * Borden discussed the Laplace- 
Poisson Mixed Equation 

(1) f{x + 1) + p(x)f(x) + q(x)f(x + 1) + r(x)f(x) = 0. 

He obtains two invariants that form a fundamental system for the equa- 
tion. These invariants are 

p(x) p(x) 

JV) = r !1 - q(.x - i). 

p{x) 

If either invariant is zero, the equation takes a simple form. For 
instance, if I(x) = 0, i.e., if 

(2) r(x) = p\x) + p(x)q(x), 
the equation (1) reduces to 

(3) f(x + 1) + p(x)f(x) = Or'^w, 

where C is an arbitrary constant. This is a linear non-homogeneous 
difference equation of the first order. Borden obtains a solution of it by 
means of the symbolic solution XG(x) of the equation 

F(x + 1) - F(x) = G(x), 

where G(x) is a known function. 

Solutions obtained in this way are purely formal, and may have no real 
significance. Borden assumes at the outset of his paper that p(x), q(x), 
and r(x) are analytic functions, but he nowhere makes an investigation to 
determine whether this hypothesis is sufficient to bestow any validity upon 
his results. 

It is the purpose of this paper to investigate in certain cases the analytic 
character of the solutions Borden obtains. To do this it is necessary to 
make use of the existence theorems for linear difference equations. We 
shall state here the one that is sufficient for several cases that arise. 

* R. F. Borden, "On the Laplace-Poisson Mixed Equation," American Journal of 
Mathematics, Vol. XLII, 1920, pp. 257-277. 
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218 Williams: The Laplace-Poisson Mixed Equation. 

Theorem A.* Consider the linear non-homogeneous difference equation 

(4) g(x + 1) + a(x)g(x) = b(x), 
where a(x) and b(x) are rational functions, so that 

a(x) = x »(a + ^+ ••• V 

bix) = x'(b<> + h+-...Y \ x \ > r. 

There exists a series g(x), of the form 

g( x ) = x-"(g + g ^+ ■X if M >0, 
-g(x) = x>(g + 9 -±+ -A, if M ^0, 

which formally satisfies (4), but ivhich in general diverges. 

There are two solutions gi(x) and gi{x) with the following properties. The 
function g\{x) is analytic save at zeros of a(x), poles of b(x), and points con- 
gruent (mod 1) to these points on the left; and it is asymptotic to the formal 
series g(x) in the right half plane. The function gi(x) is analytic save for 
poles at the poles of a(x), b(x), and the points congruent to these points on the 
right; and it is asymptotic to ~g(x) in the left half plane.] 

From the relation (2) it is possible to determine any one of the three 
functions p, q, r in terms of the other two. We shall assume that p and q 
have the form of rational functions at infinity, so that 

p(x) = x m (po + ^+ ■■■), 

(5) \ x J 

q(x) = x n (q + ^+--A, \x\>R. 

There are four cases to consider, according as n < — 1, n = — 1, 
n = 0, n > 0. 

* K. P. Williams, "The Solutions of Xon-Homogeneous Linear Difference Equations 
and their Asymptotic Form," Transactions of the American Math. Soc, Vol. XIV (1913), 
pp. 209-240. 

f In order to be assured of solutions and to know their asymptotic forms, it is sufficient 
to assume that a{x) and b{x) have the form of rational functions at infinity. In that case, 
nothing can be said of the nature or situation of the singularities of the solutions in the 
finite part of the plane. 
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§1. 
Case 1. Let n < — 1; then we have 



g—Jq(x)dx 



\ l + x~ n ~ 1+ x- n - 2+ '") 



It is seen that equation (3) is now in the form (4), with p = to, v = 0. 
It follows from Theorem A that equation (3), and therefore equation (1), 
has two solutions, analytic in general, and these solutions are asymptotic 
in the right and left half planes, respectively, to a series 

/(*) = ar-(/o+J+ •••). if to>0, 

/Or) =/„+£+ ..., if m^O, 

a; 

which can be determined by direct substitution in (3). 

The formal solution f(x) can, however, be directly obtained from (1). 
The calculation will be made only for the case m > 0. 

When we substitute in (2) the values of p(x) and q(x) given in (5), we 
obtain 



r(x) 



= x m 1 [ toj>o H h • • • ) 



the highest power of x being the (m — l)th, since the assumption n < — 1 
gives m — 1 > to + n. The equation (1) then takes the form 

f(x+l) + x m (p + ^ + ■••)/'(*) 

+ *" (q + |+ • • • )/(* + 1) + a*" 1 (rapo + -+••• W) = °- 
Assume 

/(*) = *-» (/o+^+ •••)' 
and substitute, giving 
- -^R W° + •••)-- (Po + • • -)(™fo + • • •) 



+ JL(<7o+ •••)(/o+ •••)+i(mPo+ •••)(/o+ •••) = 0, 

where we have omitted only powers of 1/x. Remembering that m > 0, 
n < — 1, we see that / , f\, / 2 , • • • can be determined step by step, the 
quantity / being arbitrary, f 

* We shall omit writing coefficients in series in 1/x where the coefficients can be deter- 
mined in terms of given quantities, and their explicit form is of no concern. 

f If m = it is found that /o is arbitrary, while some of the next coefficients, fi, fi, 
. . ., are zero, the exact number having that value depending on m and n. 
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Theokem 1. Assume the coefficients p, q, and r satisfy the following 
hypotheses; 

(1) r(x) = p'ix) + p(x)q(x). 

(2) p(x) has a pole of order m at infinity (m = 0, if p(x) is analytic at 
infinity) . 

(3) q(x) is analytic at infinity, and has a zero there of at least the second 
order. 

Then there is a formal solution f(x) of equation (1) of the for m 



/(.c) = x 



■(/.+*+•■•). 



which can be determined by direct substitution (/ is arbitrary). There are 
two solutions of (1), namely f\{x) andf 2 (x), analytic in general in the finite 
plane. Furthermore, fi(x) is asymptotic to f(x) in the right half plane, and 
f\{x) is asymptotic to f(x) in the left half plane. 

§2. . 
Case 2. Let n = — 1. In this case 



g—/q(x)dx 



\ X X 1 ) 



The equation (3) is therefore of the form (4) with v = — q * 

By application of Theorem A we deduce 

Theorem 2. Let p, q, and r satisfy hypotheses (1) and (2) of Theorem 1, 
and in addition suppose 

(3) q{x) is analytic at infinity, with a zero of the first order, and 
lim xf(x) = q . 



Then there is a formal solution f(x) of (1) of the form 
fix) = x -««-™(f 0+ f2 + ...\ 



which can be determined by substitution. There are two solutions fi(x) and 
fi{x) with properties similar to those described in Theorem 1. 

§3. 
Case 8. Let n = 0. In this case 

e -/ q (x)dx = e -gox x - qi I i _| 1_ . . 



K+-) 



* It is not necessary that v be an integer in (3). In fact the substitution 

g(x) = x'G{x) 
will reduce (4) to a normal form with v = 0, and in which the expansion of a{x) still begins 
with x*. 
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so that equation (3) is not of the form (4). It can, however, be reduced to 
that form by the substitution 

f( x ) = e-™ g ( x ). 

The new equation will be 

The application of Theorem A then gives 

Theorem 3. Let p, q, r satisfy hypotheses (1) and (2) of Theorem 1, 
and in addition assume 

(3) q(x) is analytic at infinity, with 

lim q(x) = q , lim x(q{x) — q ) = qi. 

•i'=00 .('=00 

Then there is a formal solution f(x) of (1) of the form 

f{x) = <r«*jr«i-» (f + £ + • • • V 

Solutions fi(x) andfiix) exist, and have properties similar to those described 
before. 

§4. 
Case 4- Let n > 1. We can now write 

e -fq(x)dx = f - 9 ")H. ?! (j.) ) 



where 



q w (x) = q — — - + h q n x, 

n + 1 



?«>(. 



; a .) = ar «-H^l + ?^ + _ + ...\ 



the quantities go, ?i, • • ■ having the significance given in (5). The equation 
(3) accordingly takes a form to which Theorem A is not applicable. A 
direct examination must therefore be made in order to determine whether 
a solution exists. 

Consider the associated homogeneous equation 

g(x + 1) + p(x)g(x) = 0. 

There is a formal solution 

(6) g(x) = x™*{- poe-nyx' (g + | + • • • V 

which in general diverges. There are two solutions gi(x) and g^x), analytic 
in general in the finite plane, and asymptotic to g(x), in the right and left half 
planes, respectively* 

* G. D. Birkhoff, "General Theory of Linear Difference Equations," Transactions of 
the American Math. Soc, Vol. 12 (1911), r pp. 243-284. 
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Let g(x) represent one (which one to be specified later) of these solutions, 
and put 

(7) f{x) = a,(x)g(x), 

We then have, upon substituting in (3), the equation 

(8) w(x + 1) - «(*) = C e - -i-p = <?(*). 

The discussion of this equation falls into various cases. 

Case 4a. Suppose q > 0. Let g(x) be gi(x). Then the series 

(9) a(x) = - G(x) - G{x + 1) - 

is uniformly convergent if the real part of x is positive and sufficiently large ; 
and the series is a solution of (8). The proof in all its details will not be 
given, but is based upon the following considerations.* If s is sufficiently 
large, the sih term in (9) can be written, on account of the fact that g x (x) 
is asymptotic in the right half plane to g(x), (g(x) having the form given in 
(6)), 

e -9< U (x+«-l) 

(x + *)»w-»>(- V(fi - m )*+>(x + *)•+«»+! M( - X ' 5) ' 

where M (x, s) is bounded. Now the dominant part of the numerator is, 
for large s, 



e n+l 



and the dominant part of the denominator is s m,+r+qn+1 ■ e~ m> . The dominant 
part of the expression is therefore 

e -^ L jJ" + l + ms-(«i«+r+j n+1 ) logs- 
It is seen that this is the term of a rapidly convergent series, since qo 

is positive. This is true irrespective of the sign of m, for s n+1 dominates 

s log s, for s large, if n > 0. 

Case 4b. Suppose q < 0, and let n be odd. Let g(x) be gi(x). Then 

the series 

(10) a(x) = G(x - 1) + G(x - 2) + • • • 

is uniformly convergent for the real part of x negative and sufficiently 
large; and it is a solution of (8). The proof is effected as before, making 
use of the fact that gi{x) is asymptotic to g{x) in the left half plane. 

Case 4c Suppose qo < 0, and let n be even. Neither a series of form 
(9) nor (10) will converge in this case, with g{x) representing either gi(x) 
or gi(x). This is true because s n+1 dominates s log s. 

A solution of (8) can be obtained by means of a contour integral due 

* Williams, loc. cit., § 3. 



Williams : The Laplace-Poisson Mixed Equation. 223 

to Guichard.* The integral 

. , r G(t)dt 

satisfies (8), when the path of integration passes between x — 1 and x, 
with the point x on the right, provided the path extends to infinity in such 
a way that the integral converges. 

In order to determine a choice of the path of integration we shall examine 
the nature of 

Ce-' ( "w. g W(Q = _ Ce- q(l) «>-q™{t) 
U g(t+ 1) " P(t)g(t) 

It is obvious that we can neglect </ (2) (0 and p(t), if we change the 
exponent r in the asymptotic form of g(x) to r' = r + m + q n +i- (This 
change has the effect of making q &) (t) and p{t) approach constants as t = °o .) 

Let g(t) be gi(f); then the dominating part in the denominator of 
| G{t) | is, for t large in the right half plane, 

\t mt {- p Q e~ n ) t f'\. 

If we put t = u + iv, and also t = re** this becomes 

i«i+r') logr— (m— log— p )u— m^r (t'"COS ^+r') log) — [(>«— log — j> ) cos»>+mv>sin v ] r 

Now choose the path of integration so as to make it coincident with the 

lines 

x 



<P = 



n+ 1 



at a sufficient distance from t = 0. Then t = re n+1 so that t n+1 = - r n+1 . 
It follows that the dominant part of the real portion of the polynomial g (1) (0 
will be, for t sufficiently great, the positive quantity 

ko[r" + 1 
n+ 1 

Along the contour considered it follows that 

,„,,, -,, , jrr" '— (rm C os«H-r') logi— [(»»— log— Po) cos *>+m(o sin v] T 

I G(<) | < M{r)e n+l 

where if (r) is bounded. It is seen that G(t) will approach zero more rapidly 
than e~ T as r approaches infinity, and this irrespective of the value of m, 
r', and p. 

A glance will reveal the behavior of the denominator in the Guichard 
integral on the distant part of the path of integration. Let x = £ + i-q. 

* Williams, loc. cit., § 2. 
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Then 

It is apparent that along the ray for which <p = w/(n + 1) the quantity 
just written approaches — 1, while along the ray from which <p= — ir/(n + 1) 
it becomes large. 

It follows from the considerations given that the integral converges, 
and will thus furnish a solution of the equation (8). 

We shall not inquire into the existence of other solutions or into the 
behavior of the solution in the infinite part of the plane. 

Theorem 4. Let p, q, and r satisfy conditions (1) and (2) of Theorem 1 ; 
and in addition let q(x) have a pole at infinity. Then it is possible to find a 
solution of (8), such that (7) will give a solution of equation (3), and therefore 
of equation (1). 

§ 5. 

In the identity (2) we have thus far assumed that p(.r) and q(x) were of 
given form at infinity. 

Suppose that r(x) and p{x) are given; then 

. r(x) — p'(x) 
p{x) 

If p(x) is analytic at infinity we may have any of the four cases, depending 
on the behavior of r(x). If p(x) has a pole at infinity, and r{x) is analytic, 
or has a pole of lower order than p{x), the equation reduces to Case 2, or 
possibly Case 1. If r(x) has a pole of the same order as p{x), we have 
Case 3. If p(x) has a zero at infinity and r{x) does not, the function q(x) 
will have a pole, and we have Case 4. If r(x) also has a zero at infinity 
we may have any of the four cases. 

Suppose q{x) and r(x) are given. Then 

p(x) = e - rq(z)dx [_ Ce /t(x)dx r(x)dx + C]. 

q(x) = x n >(q + ^+ ■\, 

r(x) = xn2 ( r <>+ 1 ~+ •)> \x\> R. 

In order that p(x) have the form of a rational function at infinity we 

must assume that ni < 0. If wi = — 1, we must in addition assume that 

<?o is an integer (positive or negative) such that q + m =j= — 1. In all cases 

not excluded equation (1) will come under a form already treated. 

Bloomington, Indiana, 
March, 1921. 
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